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Abstract 

We investigate imaging by spherically symmetric absolute instruments that provide perfect 

^_^ imaging in the sense of geometrical optics. We derive a number of properties of such devices, 

y^ present a general method for designing them and use this method to propose several new absolute 

• ^ instruments, in particular a lens providing a stigmatic image of an optically homogeneous region 

Q^ and having a moderate refractive index range. 

O 

•^ 1 Introduction 

^ In the last decade, perfect imaging has become one of the hot topics in optics. This was triggered by 

i i a discovery of Sir J. Pendry [1] who showed that a slab of material with negative refractive index can 

focus light to a spot much smaller than the wavelength of light. Modern metamaterials with carefully 

^ designed electric and magnetic responses |2] provided a suitable experimental background for testing 

■^ such super-resolution [3], but it has turned out that there are severe limitations due to a strong 

if~^ absorption in negatively refracting materials [3]. This started a search for other devices that could 

CN provide super-resolution but would not suffer from the disadvantages related to negative refraction. 

C It has turned out recently that such devices do exist, which was demonstrated both theoretically 

^ [5l |6] and experimentally [7] for a well-known optical device. Maxwell's fish eye [8]. Optical imaging 

T-H with super-resolution is generally called perfect imaging and the corresponding devices are called 

> perfect lenses. 

'k>( The concept of perfect imaging has become important also in geometrical optics where limitations 

\-( of optical instruments are of a different kind. There, the diffraction limit of resolution is not the 
subject of investigation, but instead one seeks to minimise or even eliminate the optical aberrations 
of the device. Remarkably, there exist devices in which this elimination can be made perfect; they are 
called absolute instruments [S]. An important subset of absolute instruments is formed by devices 
that produce images geometrically similar to the imaged objects. Within geometrical optics, such 
images are called perfect [9]. Hence the meaning of "perfect imaging" is different in geometrical and 
in wave optics. In particular, an imaging device that is perfect from the point of view of geometrical 
optics may or may not be perfect from the point of view of wave optics, and vice versa. The only 
devices that are known to produce images perfect in both senses are Pendry's slab [T] and Maxwell's 
fisheye surrounded by a mirror [HI [7] ; it remains unknown whether other lenses producing stigmatic 
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images also provide super-resolution. However, as the answer is very likely to be positive, a search 
for new absolute instruments with interesting and useful properties is highly desirable because such 
devices might provide unprecedented resolution and find applications in imaging or lithography. 

In this paper, we focus on absolute instruments and perfect imaging from the point of view 
of geometrical optics. We analyse general properties of imaging by spherically symmetric absolute 
instruments and develop a method for designing such devices. We then use this method for proposing 
several devices with interesting properties. 

The paper is organised as follows. In Sec.[2]we analyse spherically symmetric absolute instruments 
and derive several general results. In Sec. [3] we solve the inverse scattering problem for spatially 
confined rays and based on this result we develop a method for designing absolute instruments in 
Sec. |4j In Sec. [5] we give some simple examples of absolute instruments and in Sees. [6] and [7] we 
discuss devices with homogeneous regions and multiple image points, respectively. We conclude in 
Sec.il 

2 Properties of spherically symmetric absolute instruments 

In this paper we will consider isotropic spherically symmetric refractive index profiles n{r). Such 
situations have a great advantage: if some point at radius r is imaged stigmatically, then the same 
is true for all points at the same radius. At the same time, usually the object to be imaged must 
be embedded directly into the optical medium and the rays emerge in all directions from the object 
and come from all directions to the image. This is quite different from the usual imaging e.g. by a 
camera or a telescope where rays propagate more or less in one direction. 

2.1 Angular momentum and turning parameter 

It is well known (and follows for example from the analogy between geometrical optics and classical 
mechanics) that in spherically symmetrical refractive index profiles n{r) light rays propagate in a 
plane containing the centre of symmetry. This is a consequence of conservation of the quantity 
analogous to mechanical angular momentum, the magnitude of which is ^ 

L = rn{r) sin a , (1) 

where a is the angle between the tangent to the particle trajectory and the radius vector. Motion in 
spherically symmetric index distributions n{r) is therefore effectively two-dimensional, and we will 
treat it as such unless otherwise stated. 

For a spherically symmetric medium, an important role is played by the radially normalised index 
function 

N{r) = rn{r) , (2) 

which is also known as the turning parameter [TUl [H] . Consider a light ray propagating with angular 
momentum L. It follows from Eqs. ([T| and (|2| that L = N{r) sin a and hence the ray can propagate 
only in the regions where L < N{r). When it gets to the point where L = N{r), then a = 7i/2 and 
the ray propagates purely in the angular direction. Such a point is called a turning point. 

2.2 Stigmatic images and absolute instruments 

According to Principles of Optics by M. Born and E. Wolf [9], an absolute instrument is a device 
that images stigmatically a three-dimensional domain of space. A stigmatic image of a point A is a 
point B through which an infinity of rays emerging from A pass. 



We shall distinguish between two cases of a stigmatic image. In the first case, the image of a 
point A is formed at B by all rays emerging from A into some nonzero solid angle. We shall then 
call B a strong stigmatic image (or simply strong image) of A. In the second case, although there is 
an infinite number of rays going from A to B, these rays constitute just a zero solid angle. Then we 
shall call B a weak [stigmatic] image. For example, a cylindrical lens can form a whole line of weak 
images of a given point but no strong image. 

2.3 Mutual position of an object and its image 

Suppose we have an absolute instrument with refractive index n{r) that stigmatically images a point 
A with radial coordinate ta to a point B with radial coordinate tb. This means that an infinite 
number of rays emerging from A meet at B [0]. Imagine we shift the point A infinitesimally in the 
angular direction by an angle dip to a new position A' separated from A by distance r^dip. From the 
spherical symmetry it follows that the image B gets shifted by the same angle dip to a new position 
B' separated from B by r-Qdip. Now Maxwell's theorem for absolute instruments states that the 
optical length of any curve and of its image is the same, which we can apply to the lines AA' and 
BB'. Cancelling the angle dyj, we obtain 

nirx) rA = ^(^b) ^b • (3) 

Assuming that we shift the point A in the radial direction by |drA| to a point A" instead, the point 
B is shifted also in the radial direction to the point B". This follows from the fact that imaging by 
absolute instruments is conformal [U], another consequence of Maxwell's theorem: if the angle A'AA" 
is 7r/2, so must be the angle B'BB". From Maxwell's theorem applied to the lines AA" and BB" we 
then get 

n(rA)|drA| = n(rB)|drB| . (4) 

Dividing Eq. ^ by Eq. (|3]), we get 

drA , drB 



rA ^B 



with the solutions 



tb = 


= kVA 


rB = 


k 
rA 



(5) 

(6) 

(7) 

corresponding to plus and minus sign in Eq. ([s]), respectively, and an integration constant A; > 0. 

Another question is related to the mutual angular position of an object and its image. Can there 
be a situation that the object A, its strong stigmatic image B and the origin O (centre of symmetry 
of the lens) do not lie on a single straight line? Imagine such a situation. Since B is a strong image 
of A, a full solid angle of rays starting from A pass through B. We can rotate these rays around the 
line OA by some angle ip, which also moves the point B to a new position B'. The rotated rays now 
pass through B' which, due to the spherical symmetry of the lens, must also be a strong stigmatic 
image of A. This would mean that A has strong images along the whole circle that is obtained by 
rotation of the point B around the axis OA, which is clearly impossible. So we must conclude that 
the points A, B and O lie on a single straight line. This way the image B is either on the exactly 
same side from O as is A or on the exactly opposite side. 

Note that this argument is not valid in two dimensions where the mutual position of a point and 
its image in a rotationally-symmetric absolute instrument is less restricted. For the same reason this 
argument is not valid in 3D for weak images. 



2.4 Value of the constant k in Eq. ([6 



An interesting question is related to the possible values of the constant k in Eq. ^. We conjecture 
that the only possibility is A; = 1. Although we have not been able to show this in the completely 
general situation, the proof can be given in two practically important cases. For the first case, assume 
that the constant k is the same for all the points that are imaged by the lens. But if B is an image 
of A, then also A is an image of B, so r^ = kr^ must also hold along with Eq. ([6]), from which then 
k = 1 follows. 

The assumption of the second case is that among the rays that contribute to the image B of a 
point A there is also the ray for which A is a turning point. In other words, we assume that a ray 
starting from the point A in the angular direction makes it to the image of A formed at the point B. 
Since Eq. ^ can be rewritten as N{r\) = A^(rB), it follows that the point B is a turning point for 
the same ray as well. Next we use Eq. (Is]) with the plus sign from which Eq. ^ has been derived, 
i.e., drx/rA = dre/rB. From this equation combined with N{rA) = N{rB) then follows 

ta —. = re —, . 8 

drA dre 

We see that the derivative dN/dr has the same sign at both r^ and r-^. However, if A; 7^ 1, this is a 
contradiction with the fact that the same ray has its turning points at A as well as B. To see this, 
suppose that k > 1 {in the opposite case we can relabel the points A and B), so r^ > r^- Then if 
the derivative dN/dr is negative at r = ta, then for the ray with angular momentum L = N{rA) the 
region r > ta is inaccessible because in this region L > N{r), which is impossible, and so the ray 
cannot make it to B. On the other hand, if the derivative dN/dr is positive at r = ta, then it is also 
positive at r = r-g. For a similar reason then the ray with angular momentum L = N{rA) = N{r-B) 
cannot make it to A, which is a contradiction. Hence the only possibility is ta = tb which then 
implies k = 1. 

We thus see that there are just two cases of imaging by spherically symmetric absolute instru- 
ments: either the image is given by spherical inversion (since ta^b = k) of the object, possibly 
combined with some rotation, or it is congruent with it. We have arrived at a slightly stronger 
statement than derived in [H] for absolute instruments in general. We will see in Sec. [ijthat Eq. ([T]) 
corresponds to the generalised Maxwell's fish eye; all other spherically symmetric absolute instru- 
ments correspond to Eq. ^ with k = 1, so they give images congruent with the object and their 
magnification is unity. This means that the imaging by such a device is perfect in the sense of 
geometrical optics [9]. 

3 Inverse scattering problem for spatially confined rays 

In the inverse scattering problem in mechanics, the task is to determine an unknown potential from 
the scattering angle which is a known (e.g. measured) function of the impact parameter. The 
problem was solved for central potentials in 1953 by O. B. Firsov [T2] and the analogous problem in 
optics was solved by K. Luneburg [13]. In the situations considered there, the particle or light ray 
incides from infinity, undergoes scattering and leaves for infinity again, so the motion is not spatially 
confined. However, the situation where the motion is restricted to a finite region of space may be 
equally important for design of absolute instruments and perfect lenses. The inverse problem can 
be formulated and solved also in this case; the solution for mechanical motion was given without 
derivation by V. N. Ostrovsky in [Hj. In the following we derive the inversion formula for confined 
motion in the optical case. 



To derive the inversion formula, we will assume that the function N : r \-^ N(r) in Eq. ^ 
is increasing for r < tq and decreasing for r > tq with some radius r^ > 0. (There are obviously 
more general profiles N{r) and the inverse problem can be solved for them as well, but we will not 
consider them here.) Then the function N{r) has a global maximum Lq = maxA^ = A^(ro) at the 
point r = tq. The inverse function r : N \-^ r{N) is multi- valued, i.e., has two branches r±. One 
branch r_ : N \-^ r_{N) maps into the inner region [0,ro], and the other branch r^ : N \-^ f'+i^) 
maps into the outer region [ro,oo). In this case we can invert the inequality L<N{r) into a double 
inequality 

r_(L) <r <r+(L), (9) 

which explicitly specifies the allowed confined region. In other words, for angular momentum L < Lq, 

there are two turning points r±{L). Moreover, a light ray with angular momentum L = Lq will 

propagate along a circular trajectory with radius tq, while angular momentum L > Lq is forbidden. 

For the subsequent calculations, it turns out to be more convenient to work in terms of a new 

coordinate x = Inr rather than the radius r. We will also introduce corresponding notation x± = 

\nr±, xq = Inro, etc., in the obvious fashion. At a general point of a ray trajectory, the derivative of 

the polar angle ip is 

dip dip L 

dx 



dr 



tana 



^/m{r)-L^ 



(10) 



where we have used that L = N{r) sin a. With the help of Eq. (10), the increment of the polar angle 



corresponding to motion between r_(L) and r+(L) can be written as 

'r-(L) r ^ 



Av^(L) 



■y/A^ 



L2 



dx 



(L) ,/N^W^L^ 
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We shall call Aip the turning angle; it expresses the change of the ray direction between two conse- 
quent turning points. 

The task of the inverse scattering problem is to find the refractive index n{r) [or equivalently, the 
function N{x)] from the known turning angle A(p{L). It can be solved in a similar way to finding the 
ID potential from the known period of oscillations as a function of energy [IS]. We divide A(f{L') 
by -\/L'2 — L^, where L' is an integration parameter, and integrate with respect to L' from L to Lq: 



Lo 



AipjL') dV {TT| 



Lo 



x+(L) 

x.{L) 
x+(L) 

x_(L) 
x+(L) 

x^{L) 



-{L') 



dx 



L'dL' 



N(x) j^, ^^/ 



L'^^W- 
dx 



L2 



arcsm 




(12) 



where we have inverted the order of integration and changed the integration limits appropriately, see 
Fig. [I] Equation ( 12 ) yields the inversion formula 

r^(L)_2 r^o A^{L')dL' 



In 



which solves the inverse scattering problem. The mechanical equivalent of Eq. (13) was presented in 
without derivation. 
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Figure 1: The change of integration hmits in Eq. (12) illustrating that J^° J^, ,j^, 
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4 Method for designing absolute instruments 

In the following we will describe a general method for designing absolute instruments. This method 
was sketched in p^ but was not given explicitly. 

Consider the situation when the turning angle Aip is independent of L and equal to tt/tti, m G M. 
This is a practically important case; if, for instance, m = p/q with p, q coprimes and p even, then 
the trajectory will be symmetric with respect to rotation by vr around the origin in the plane of 
propagation. A strong stigmatic image of a point A at the position r will therefore be formed at 
— r. If m = p/q with p odd, a strong stigmatic image of a point A will be formed at A itself; the ray 
arrives there after encircling the origin q times. This will be shown in Fig. |5| 

For Aip{L) = Tr/m we get from Eq. (13) ln[r+(L)/r_(L)] = (2/m) arcosh (Lq/L), which can be 
expressed as 



L 



r-{L) 



m/2 



+ 






m/2 



(14) 



Now comes the key step of our derivation. We define a function /(r) such that 



f{r) 



r+{N{r)) 



for r < tq 
for r > Tq 



(15) 



The function /(r) is hence defined such that for a given lower turning point r_ it produces the upper 
turning point r+ corresponding to the same angular momentum and vice versa: 



r± = f{r^] 



(16) 



We can say that the point r^ is dual to the point r_ and vice versa. It also follows from this definition 
that applying the function / twice yields the original value, /(/(r)) = r, and therefore the graph of 
/ is symmetric with respect to the axis of the first quadrant. The graph intersects this axis at the 
point r^ = r_ = rQ, which corresponds to the circular ray trajectory. 

With the function /(r) defined this way, we can express L from Eq. (14) and either keep r_ and 
replace r_|_ by f{r-)^ or keep r+ and replace r_ by f{r^). Taking then advantage of the fact that 
L = A^ at r-t and omitting the index of r±, we get 



N{r) = 2Lc 



f{r) 



m/2 



+ 



f{r) 



m/2 



(17) 



Finally we express n{r) 



N{r)/r as 
n{r) 



2Lo 



fir) 



m/2 



+ 



m 

r 



m/2 



:i8) 



Equation (18) provides a powerful tool for designing absolute instruments. For any chosen func- 



tion /(r) satisfying the condition above and a suitable value of m it gives a refractive index profile 
with focusing properties, i.e., an absolute instrument. We will demonstrate this on several known 
examples first and then proceed to new devices. 



5 Examples of absolute instruments 

In all the following examples we assign unit radius and unit refractive index to the circular ray, so 
rQ = Lq = 1. In the first example the object and its image are related by spherical inversion, in the 
other ones they are congruent, as discussed in Sec. |2} 

• Generalised Maxwell's fish eye 

Consider an absolute instrument in which the position of object and image is given by Eq. (JTl). 
Apparently, if ta is a turning point, then tb must also be a turning point, which leads to 



r_|_r_ 
f{r) 



1 (the case of a general k can be obtained easily by spatial scaling) and hence 



k 



1/r. Then Eq. (18) yields 



n{r) 



-l—m 



_|_ J, 



m+1 



(19) 
Light rays in this lens are 

(20) 



which is the generalised Maxwell's fish eye profile discussed in 
shown in Fig. M (a) for m = 1/2. 

For m = 1 we get 

^^^) = T^^ ' 

which is the well-known Maxwell's fish eye refractive index profile [8]. The trajectories are 
circles intersecting the unit circle at two opposite points. Light rays are shown in Fig. |2] (b). 

Lunehurg lens profile 
Take m = 2 and /(r) 



V2 - r2. Then Eq. (18) yields 



n{r) = V2 - r2 . (21) 

For r < 1 this coincides with the well-known refractive index of Luneburg lens [13] which has. 



however, n = 1 for r > 1. The index (21) corresponds to Hooke potential in mechanics and ray 



trajectories are ellipses centred at the origin pjj, see Fig. |3] (a 

Eaton/Minano lens profile 

Take m = 1 and /(r) = 2 — r. Then Eq. (18) yields 



n 



(r) = v/(2/r) - 1 . 



(22) 



This is the well-known index profile of Eaton or Minano lens ^7\ [18] which have, however, 
n = 1 for r > 1 or r < 1, respectively. The index ( [22] ) corresponds to elliptic motion in Newton 
potential in mechanics and ray trajectories are confocal ellipses with focus at the origin and 
with the main semiaxes of unit length [11], see Fig. [3] (b). 
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Figure 2: Ray trajectories in (a) generalised Maxwell's fish eye with m = 1/2 and (b) Maxwell's fish 
eye. The source and its image are denoted by A and B, respectively. Here as well as in all subsequent 
figures, the black circle corresponds to the circular trajectory at r = ro with the maximum possible 
angular momentum Lq. 



Maxwell's fish eye mirror 

Take m = 2 and /(r) = 1 for r < 1 = rg. From the symmetry of the function / it follows that 
it is undefined for r > 1; the rays cannot get beyond the unit circle, which means they must 



be reflected there. Substituting into Eq. (18), we get for r G [0, 1) the index given by Eq. (20) 



This is the so-called Maxwell's fish eye mirror discussed in detail in |5] with rays shown in 
Fig.i 

Lenses generated by a sample function f{r) 

To show the generality of our method, we choose some arbitrary function /(r) with the restric- 
tion that it is symmetric with respect to the axis of the first quadrant. Let us choose 



fir) 




for r < 1 
for r > 1 



(23) 



see its graph in Fig. p^ (a). Using different values of m, we get different absolute instruments. 
The rays in some of them are shown in Fig. M (b) - (d) for m = 2, m = 5/2 and m = 2/3. 



5.1 Lenses with mirrors and index discontinuities 

On the example of Maxwell's fish eye mirror we have seen that an interval of r on which the function 



/(r) from Eq. (16) was constant corresponded to a spherical mirror. As we will see, this is quite a 
general feature. 

Consider the situation when the function /(r) is constant and equal to c > ro on some interval 
(a, 6) with b < Tq. In other words, the upper turning point is the same for different lower turning 
points and hence for different angular momenta; the corresponding rays cannot get beyond r = c, 
which means there occurs reflection at this point. Now there are two cases to be distinguished. In 
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(b) 



Figure 3: Ray trajectories in the profile given by (a) Eq. (21) corresponding to Luneburg profile 



and (b) Eq. (22) corresponding to Eaton/Mifiano profile. The optical medium is shown in light blue 
colour; the refractive index goes to zero at the edge of the medium. In Eaton/Minano profile the 
image B coincides with the source A. 



the first case f{r) = c holds also for < r < a, which means that a reflection occurs for all rays 
that reach the point r = c; this corresponds to a perfect mirror placed at c. In the second case, the 
value /(r) gets larger than c for some r < a. This means that rays with small angular momenta can 
penetrate beyond c; the reflection at c mentioned above must then be total internal reflection caused 
by a refractive index discontinuity at r = c. That such a discontinuity indeed exists follows from the 
fact that if /(r) = c on the interval (a, b), then /(r) must be discontinuous at r = c and so must be 
n{r) according to Eq. (18). 



If the interval of constant / lies above tq then there will be a total reflection on a sphere "from 
the outside" . This is illustrated in Fig. [6] where there is a mirror reflecting perfectly from the inside 
and a jump of refractive index reflecting totally some rays from the outside. 

Spherical mirrors can be used with a great advantage to reduce the size of the lens and also 
the range of the refractive index. An example is Maxwell's fish eye mirror discussed above that has 
equally good focusing properties as the original Maxwell's fish eye but its size is reduced from infinity 
to a unit disc (or sphere) and the refractive index range is reduced significantly from the interval 
(0,2] to just [1,2]. 



6 Absolute instruments for homogeneous regions 

In many absolute instruments such as Maxwell's fish eye, the optical medium fills the whole space 
and the object to be imaged must therefore be inserted into such an inhomogeneous medium. It is 
desirable to find optical devices that provide images of optically homogeneous spatial regions, i.e., 
regions with a uniform refractive index. Even if the refractive index of such a region differs from 
unity, this is still an advantage because one can fill this region with a suitable liquid and place the 
object in it. 

Interestingly, until recently the only known devices providing stigmatic images of optically homo- 
geneous 3D regions were plane mirrors or their combinations [5]. This has changed by an excellent 





(b) 



Figure 4: Ray trajectories in (a) Maxwell's fish eye mirror an (b) Minano lens. Optically homogeneous 
region is shown in white colour. 



work of J. C Miiiano [18] who noticed that some well-known optical devices such as Eaton or Luneb- 
urg lens [T^ are in fact absolute instruments, providing stigmatic virtual images. An even more 
important result of [IB], however, is a a device now called Minano lens that provides real images of 
homogeneous 3D regions. 



6.1 Minano lens 

In Miiiano lens the homogeneous region is the unit disc with unit refractive index and the turning 
angle is Ays = 7r/2, which corresponds to m = 2. To derive the refractive index outside the unit disc 



by our method, we will again use Eq. (18), but first we have to determine the function /. The lower 
turning point for a given L is simply r^ = L because n = 1 in this region, the upper turning point is 
r+ = f{r_). Substituting this into Eq. (14), we get a quadratic equation p{r_) — 2/(r_) + r^ = 
for f{r_) with the solution 

-r2 (24) 



fir- 



1 + 



(we have taken the larger root so that r+ > r_). Inverting this expression to get the function / also 

\/2rj^ — r\. Now we can combine Eqs. (18) and (24| with Lq = 1 



for r I , we find r_ 



/(r+) 



to find the refractive index outside the unit disc. This gives precisely the expression (22), so the 
refractive index of Minano lens is 



n(r) 



1 



vW^ 



for r < 1 
for r > 1 



(25) 



and the ray trajectories are shown in Fig. |4] (b) 



6.2 Modified Minano lens 

Despite its elegance and nice properties, Mifiano lens has a disadvantage: its refractive index ranges 
from unity all the way to zero at r = 2, which is very difficult to realise practically. Since any 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

(a) 



(b) 





(d) 



Figure 5: Graph of a sample function /(r) according to Eq. (23 ) and the corresponding ray trajectories 
for (b) m = 2, (c) m = 5/2 and (d) m = 2/3. In (b) and (d), there is a strong image of A at B while 
in (c) the images B, C, D, E of A are just weak. 



refractive index profile can be multiplied by a real number without affecting the lens performance, 
we can define a measure of the index range as the ratio of its largest and smallest value 



V 



n 



/rir 



max/ "-niin 



(26) 



For Minano lens rj = oo. It would be very desirable to modify the lens somehow to make t] finite. 
y, this is possible, as we will show now. Imagine we use the function / for r < 1 according 



Fortunate^ 



to Eq. (24), but on the interval [0, b] with < 6 < 1 we replace it by a constant value c = l + -\/l — b"^. 



As we have seen, this corresponds to placing a mirror at r = c. The function / this way becomes 



fir) 



see Fig. [7] (c), and refractive index 



n{r) 



1 + v^ 
V2 






'2c/(r2 + c2) 



2/r- 1 



for r < b 
for 6 < r < 1 
for 1 < r < c 



for r < b 
for 6 < r < 1 
for 1 < r < c 



(27) 



(2^ 
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0.0 0.5 



(d) 



Figure 6: (a) Ray trajectories and (b) imaging by an absolute instrument with a spherical mirror at 
r = 2 and a jump of refractive index at radius r = 1. There are two types of ray trajectories differing 
by whether they do or do not penetrate inside the unit circle, (c) The corresponding function /(r) 
and (d) refractive index. 



see Fig. u\ (d). The largest and smallest values of the index occur at r = and r = c, respectively, 
and yield the ratio r] = 2/6. Choosing b not too small, one can get a very moderate index range of 
the lens; the price to pay is that the size of the homogeneous region is reduced. Rays in modified 
Miiiano lens and its imaging properties are shown in Fig. [T] (a) - (b). 



6.3 A lens for designing a magnifying absolute instrument 

Recently we proposed another type of lens for imaging homogeneous spatial regions [19] that has a 
finite ratio of the largest and smallest refractive index. The homogeneous region is a unit sphere 
and there is a spherical mirror at radius R > 1. The refractive index between the two spheres is 
chosen such that a ray emerging in the direction of a unit vector u from a point A located at f^ in 
the homogeneous region incides on the mirror at the point Ru. This ensures that the ray after the 
reflection from the mirror passes through the point B at re = — ^A where an image of A is formed, 
see Fig. |8| Another property of the device is that, apparently, all mutually parallel rays propagating 
in the homogeneous region are focused to a single point at the mirror. 
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b) 



n(r) 



2/c 



b/c 




c) 



b 1 

d) 



c 2 



Figure 7: (a) Ray trajectories and (b) imaging by the modified Miiiano lens. In contrast to Mifiano 
lens, there is a spherical mirror at radius c and an inner inhomogeneous region of radius h. (c) The 
corresponding function /(r) and (d) refractive index n{r). The dashed parts of the graphs on the 
intervals [0, h] and [c, 2] correspond to the original Miiiano lens. 
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Figure 8: Rays in the lens imaging a circular homogeneous region employing a spherical mirror. The 
homogeneous region is encircled by a dotted line. Unfortunately, the refractive index performing 
such an imaging perfectly does not exist as we show in Appendix. 

The refractive index calculated for this lens was based on numerical solution of a specific integral 
equation that does not have an analytic solution. Unfortunately, it has turned out later that the 
equation does not have solution at all, so the proposed lens is in fact not an absolute instrument. The 
rays emerging from A hence do not pass exactly through B but there is some unavoidable aberration. 
However, numerical simulations reveal that this aberration can be very small if R is not too small. 
For this reason the lens could have a practical importance, and we are therefore mentioning it here 
although it is not really an absolute instrument. The proof of non-existence of the solution of the 
integral equation is given in the Appendix. 

7 Lenses with multiple image points 



So far, we have considered just situations when the turning angle A(p{L) in Eq. (11 ) is constant for all 
angular momenta L from zero to Lq. However, another class of interesting absolute instruments can 
be obtained by taking Aip{L) piecewise constant on different intervals of L. In this situation, however. 



the refractive index cannot be found analytically with the help of Eq. ( 13 ) but must be calculated 



numerically. The simplest situation corresponds to two values of turning angle: Aip{L) = ir/rrii for 
L E [0,Li] and Alp(L) = n/m2 for L E [Li,Lo]. Fig. |9] (a) - (b) shows an example of this lens for 
mi = 1 and m2 = 2. Rays emerging from the point A with angular momentum larger than Li meet 
first at the point B where they form an image and then continue back to A. On the other hand, rays 
with L < Li go back to A directly without passing through B. This way, the image B is formed just 
by some rays while the image A is formed by all rays. 

Since there is not a smooth transition from the first class of rays (those with L > Li) to the 
second class {L < Li), the optical path from A to B does not have to be the same for rays of the 
two types. In other words, the principle of equal optical path [9J does not apply to this type of lens. 

Another example of such a "bifocal lens" is shown in Fig. M (c) - (d) for mi = 2 and m2 = 4. 
In the two-dimensional version of this problem, some rays emerging from A form an image at C or 
D and then continue to the image B. Other rays go to B directly without passing through C or D. 
If the lens is three-dimensional, a strong image will be just at B (and then at A again), but weak 
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images of A will be formed on the circle made by rotation of the points C and D about the axis AB. 
For coherent light waves emitted from the source at A, the inequality of the optical paths from A to 
B for the two ray types could lead to interesting diffraction patterns around the image point B if the 
two corresponding waves interfere destructively. 





b) 





c) 



d) 



Figure 9: Ray trajectories in "bifocal" lens with Lq = I, Li = 3/4, /(r) = \/2a'^ — r^ and (a) - 
(b) nil = 1, ""^2 = 2 and (c) - (d) mi = 2, m2 = 4. The dashed circles mark the radii of turning 
points corresponding to angular momentum Li. Rays confined to the ring between these circles have 
turning angle Aip = TT/m2 while rays that get outside the ring have Aip = tt/itli. 



8 Conclusion 

In this paper we have derived a number of results for imaging by spherically symmetric absolute 
instruments. We have shown that images created by such instruments are either congruent with the 
object or related to it by spherical inversion. We have also proved that the mutual position of a 
point and its strong image is quite restricted in 3D: the straight line connecting the two points must 
intersect the centre of the lens. 
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Further, we have developed a general method for designing absolute instruments via solving the 
inverse problem for finite motion. We have shown how this method can be used for designing already 
known as well as a number of new absolute instruments. The modified Minano lens we have proposed 
is particularly appealing because it produces a real stigmatic image of a homogeneous region with 
a reasonable range of refractive index, and we believe that this lens can find important practical 
applications. Whether the imaging devices discussed in this paper can provide super-resolution 
remains an open question, but in our opinion the answer is positive. 
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Appendix: No-Go Theorem for Spherically Symmetric Lenses 



In the following we show that the refractive index with the appealing properties described in Sec. 6^ 
does not exist. This gives a negative answer to the following question: is there a spherically symmetric 
medium that focuses parallel rays in a spherical cavity via a shell-shaped lens into a single point? If 
the answer were positive, this would constitute a new interesting and practically important absolute 
instrument as well as a useful building block for other optical devices [IH]. Although the actual 
answer is negative by itself, it still useful and likely to influence future developments in the subject 

m- 

Consider parallel rays located in an inner spherical region r<r^ with a constant refractive index 
n^ = — >0 that are focused by the medium to a point P situated at some finite distance R>r^: from 
the centre O. We will try to determine an appropriate spherically symmetric profile N{r) = rn{r) in 
the intermediate shell region 

n < r < R (29) 

with inner and outer radius r^, > and R < oo, respectively. We will assume that all rays are 
outgoing, i.e., the radial position of a light ray r(t) is monotonically increasing as a function of time 
t, so no turning points are allowed. In particular, the outer region r>R will play no role, and we 
might as well assume that the focal point P lies on the outer rim of the device. Suppose that the 
parallel rays move horizontally to the left and are focused to the focal point P with polar coordinates 
(i?, vr). As before, it is convenient to introduce the logarithm x = lnr of the radial coordinate r. 
Similarly, we shall write a;* = lnr*, X = \nR etc. in an obvious notation. 

Clearly A^^, is the upper bound for the angular momentum L inside the inner region r < r^:. 
However, we shall only require that the focusing device actually works for rays with L e [0,Lm), 
where Lm G (0, A^*]. A ray with angular momentum L emerges horizontally to the left from the 
inner rim at a polar angle n — arcsin -^ and when it arrives at P, the polar angle becomes vr, so 
the corresponding change of polar angle is arcsin^. On the other hand, the same change can be 



calculated by integrating Eq. ( 10 ) from r* to R, so we arrive at the following integral equation: 



VL e [0, Lm) : 7 arcsin # = f^ , '^' = f ^ , ^^- . (30) 

We next redefine (normalise) the quantities N' := -^, L' := -^, L'^ := ^, etc., by dividing with the 
constant A^*. Suppressing the primes from now on, we arrive at the main integral equation pTOj, 
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which serves as our starting point. The solution function N = N{r) is required to be bounded 

3N^ Vr e [n, R]: < L^ < N{r) < iV„ < oo (32) 



for physical reasons, and Lebesgue measurable [21] in order for the integral (31) to be defined. Since 



the integrand (31) is non-negative, the integral is always defined, although perhaps infinite. The goal 



Note that the fc'th moment f 



bounds (29) and (32). 



R dr 1 

r N''{r) 



is now to investigate the integral equation (31 ), and show that it has no bounded solutions N = N{r). 



is well-defined and finite for any integer k because of the 



Lemma 8.1 (Reformulation in terms of odd moments) The integral equation (31) is equiva- 
lent to 

^ ' (33) 



X 



R dr 1 



^^ for nGNo = {0,l,2,3,...}. 



Ir, r Af2"+i 2n^ 

Proof: The Taylor expansion of the arcsin(L) function at L = is 



arcsm(L) 



Eoo 1 

n=0 2n+l 



2n 



n 



L\2n 



(f) 



(34) 



On the right-hand side of Eq. (pl|, a binomial expansion produces 

r^ dx _ fX dx A _ ( I±\2\~2 _ f^ dx ■spoo / 

Jx, ^N'^-L'^ ~ Jx, N \^ VAf/ / ~ Jx, N 2^n=0 1 

_ fX dxsr^oo ( \ ( Jl_\'^"' — Y^°° I 

~ Jx, N 2^n=0 \ fi j \2n) ~ 2^n=0 I 



n 
2n 
n 



fX dx / L \2" 
Jxt N \2n) 



(35) 



In the last equality of Eq. (|35|, we were allowed to exchange integration and summation order. 



because each term is positive or zero (Tonelli's Theorem). In the third equality of Eq. (35), we have 
used the identity 

' 2^ ^ {-\r . (36) 



_1 

2 

n 



n 



Comparing the Taylor-coefficients on the left-hand side rt34| and right-hand side (35) of Eq. (31) 



yields the reformulation ( 33 ) . D 



VLgC 



(Reformulation in terms of new L) The integral equation (31) is 
■ e'' = Xf^(l + 2(^)^)exp[(^)^] ^ (l + £)r^exp[(f)^] 



a=l 



(37) 



Proof: 

calculate 



The integrals in Eq. (37) are well-defined due to the bounds (29) and (32). We 



, 2n 



te3 



r r^ dx \(I^\2] _ J r^dxv^oo }_(I^y^ _ Y^oo 

^ Jx, N ^^^^ Wn' J ^ Jx, N Z^n=0 n! \n) Z^n=0 n\ 

„oo 1 L^"+l _ 'Spoo J_ cL -tp p2n _ r^AP\^°° J_/72n 
l^n=0 n\ 2n+l Z^n=0 n! JO ^^ ^ JO ^^ Z^n=0 n!^ 



hixX^r^' 



(3^ 



In the second equality of Eq. (38), we use Tonelli's and Fubini's Theorems [21] with x-independent 
majorant -M- exp (-M-)^ to justify exchange of integration and summation order. Tonelli's and 



Fubini's Theorems are also used in the fifth equality of Eq. (38) with majorant e' ' . Now differentiate 



both sides of Eq. (38) with respect to L to obtain Eq. (37). Note that Eq. (38) is trivially satisfied for 
L = 0, so we do not lose information when differentiating. Therefore we can also run the argument 
backwards. D 



Thus we have three equivalent conditions, Eqs. (31), (33) and (37), that a solution N = N(r) 
should satisfy. 
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Lemma 8.3 // there exists a solution N = N{r) (not necessarily monotonically increasing as a 
function of the radius r), then there also exists a monotonically increasing solution N = N{r). 



Sketched proof: Lemma 8^ is basically the observation that the integration variable 

X = Inr only enters implicitly via the function N = N{x), and that the Lebesgue measure dx is 
translation invariant. D 



From now on we can and we will make the following assumption |8.4| without loss of generality. 
Assumption 8.4 The solution N = N{r) is a monotonically increasing function ofr. 



At this point, we introduce a technical assumption 8^ in order to proceed. 

Assumption 8.5 The inverse solution r = r{N) exists and is differentiahle with Lebesgue measurable 
derivative. 



The Assumption |8.5| implies that one can define a Lebesgue measurable density 



(39) 



Let us call the definition domain of the inverse solution r = r{N) for [A*, A^]. In other words, 
r(A,) = r, and r(A„) = R. 



Lemma 8.6 (Reformulation in terms of test functions) Under the assumptions 8.4 8.5, the 
integral equation (31) becomes equivalent to 



WV e C7-((0,oo)) : /;;XA)dA ^ = -r/(l) . 



(40) 



Remark: Here C^((0, 00)) denotes the set of infinitely often differentiahle functions rj defined 
on the open interval (0, 00), and such that 77 has compact support in (0, 00). Compact support means 
that the function rj is assumed to vanish identically in whole neighbourhoods around of A = and 
A = 00. It is therefore natural to extend r] smoothly to the closed interval [0, 00] by assigning to 77 
the values 77(0) = = 77(00) at the end points A = and N = oc. 



Proof: When we substitute the inverse solution r = r(A), equation (37) becomes 



(i+£)/;:pw^-p[(f)^] 



a=l 



(41) 



Next perform the elementary substitution u = 1/A with limits z/^ = l/N^, and z/^, = 1/A^,. Fur- 
thermore, multiply both sides with e~^^^^ , where /x > is a positive parameter. Then 



,L2(l-/.2) 



(1 + £) /:>(-)d- e 



L2((q,^)2_^2) 



a=l 



(42) 



Recall that the Dirac delta distribution S{x) has the Fourier integral representation 5{x) = J^ 2^ ^'^^■ 
By integrating L^ = ip along the imaginary axis in Eq. (42), one gets 



^(1-/^') = (i + £)/::p(^)d^^(M^-^^; 



a=l 



By using the Jacobian formula for the Dirac delta distribution 



sifix)) = E 



1 



f{xo) = 



6{x — xq) 



(43) 



(44) 



and multiplying both sides with 2, one gets 



«=i 



(45) 



where we have assumed that ;U > is positive. Thus for a test function t] G C^((0, oo)), one 
calculates 



Vil) = /,°^d^ r^(^)<5(/i - 1) ^ {l + £)j2piu)duJ^-df,'-f6{au-^^) 



(1 + £) /:>(-)d^ 



Ti{au) 






a=l 



^uyi^)^^ 



ri{ciu) 
au 






a,=l 
a=l 



(46) 



Now translate (46) back to the A^ = l/u variable to obtain Eq. (40). D 



Lemma 8.7 Eq. (40) has no solutions for p that respects the bounds (32) on N. 



Proof: The Fundamental Lemma of calculus of variation (in the strengthen version of du Bois- 
Reymond) [22] shows that p must be a constant up to contributions that vanish almost everywhere. 
(In particular, we stress that it is not enough for p to be only piecewise constant.) Thus one may 



pull the density p outside of the integral (40), and integrate to get 



\/veC^iiO,oo)): p(r/(iVj-r^(iV,)) 



-^(1) • 



(47) 



Collapsing limits A^* = A^^ are clearly not a solution. Assuming A^^, < A^^, equation (47) has two 



solutions, (p = l, A^=K = 1, A^m = oo), or (p = — 1, A^^, = 0, A^m, = l)- However, none of these two solutions 



respect the bounds (32) on A^, and the latter is actually monotonically decreasing. D 
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